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ABSTRACT 

We propose a perturbative improvement of the hierarchical approximation for 
gaussian models. The procedure is based on a relabeling of the momenta which 
allows one to express the symmetries of the hierarchical model using a simple 
multiplication group. The representations of this group are used to expand the 
action. The perturbative expansion is treated as a problem of symmetry breaking 
using Ward identity techniques. 



1. Introduction 



In many respects, gauge theories provide a satisfactory framework to describe 
the interactions among elementary particles at energies accessible with existing 
colliders. Nevertheless, the perturbative methods which are spectacularly accurate 
for QED at low energy are not adequate to take into account large distance effects 
among strongly interacting particles. The renormalization group (RG) method is 
an essential tool to understand and describe these phenomena and, more generally, 
the infra-red behavior of a field theory. However, the practical implementation 
of this method requires approximations which are usually difficult to control or 
improve. The goal of this paper is to discuss the improvement of such an approxi- 
mation in a very simple scalar theory. 

In his remarkable analysis of the partition function of the scalar field theory, 
K. Wilson '^'has shown that by replacing some of the contributions by order of 
magnitude estimations, it was possible to obtain a very simple renormalization 
group transformation called the approximate recursion formula. This recursion 
formula played an important role in the development of the RG method because 
the basic ideas (fixed points, relevant directions..) were not too difficult to work 
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out in detail ' in this simplified RG transformation. Similar recursion formulas 
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hold exactly for hierarchical models. For this reason, we call the approximation 
made by Wilson the hierarchical approximation. 

r^i 

An interesting feature of the hierarchical models is their large group of sym- 
metry. An appropriate use of these symmetries allows, in the case of Ising spins, 
to cut down the time necessary to calculate numerically'*' the partition function 
with 2"' sites, from T — 2^" to (log(T))'^, without making any approximations. In 
our attempt to improve the hierarchical approximation, we would like to keep the 
largest possible subgroup of unbroken symmetries at each stage of the calculation. 

The improvement of the hierarchical approximation is in principle straightfor- 
ward: we have to restore the details erased in Wilson's discussion!^' However there 
are several practical problems which will be encountered. First, there is a book- 



2 



keeping problem: we need to write all the corrections. Second, the is a priority 
problem: we would like to know which corrections are the most important in order 
to calculate them first. Finally, there is the feasibility of the calculation itself. 

We present here a method which overcomes these problems in a natural way. 
For simplicity, we have restricted the discussion to a gaussian model on a finite 
one-dimensional lattice with 2*^ sites. The Fourier modes of the scalar field are 
denoted where = and k are integers used to express the momenta in 
1^ units. These integers are understood modulo 2^'' in the following (periodicity in 
momentum space). The action reads 



where g{k) is even, real and positive . 

We proceed in three steps. First, we relabel the momenta in a way which 
is convenient to read their "shell" assignment in Wilson's discussion (section 2). 
This relabeling allows us to introduce a group of transformation whose orbits are 
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precisely these shells. The representations of this group are known and can be 
used to expand the kinetic term , i.e., the function g{k) in Eq. (1), for each of the 
shells. This solves the bookkeeping problem. Nicely enough, the classification of 
the representations of the group mentioned above comes with an index indicating 
its resolution power (called degree of ramification). This naturally provides the 
successive orders of our perturbative expansion. We then show that if we only 
retain the trivial representation in the expansion, we obtain the hierarchical ap- 
proximation (section 3). In this limit, the group of transformation is a symmetry 
of the action which can be identified with the symmetry group of the hierarchical 
model mentioned above. The improvement of the the hierarchical approximation 
can thus be treated as a symmetry breaking problem. We describe the first cor- 
rections using the familiar apparatus of the Ward identities in section 4. 

We emphasize that this procedure takes advantage of the symmetry present at 
order zero in an optimal way. We have checked in simple examples that the largest 
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k=l 
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corrections correspond to the smallest breaking of the symmetry. The extension 
of this construction for D-dimensional models is rather straightforward and men- 
tioned at the end of the paper. The use of this scheme for interactive theories 
( Landau- Ginzburg or Ising type) is presently under study and briefly discussed at 
the end of section 3. 



2. A Convenient Relabeling of the Momenta 

In this section, we introduce a relabeling of the momenta k appearing in Eq. 
(1). For this purpose, we use a set orthonormal functions which is a discrete version 
of the Walsh system ''''motivated by Wilson shell decomposition (see section II of 
Ref.[l] for detail). We first define 



uii k = 2'^-2 
ou* iik = -2'^-2 
. otherwise 



(2) 



and 



*i(A;) = 1 - ^o{k) 



(3) 



with the notation u = 



i+i 

2 ■ 
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For a given integer a = ao + ai2^ + 022^ + + an-i2^ ^ with a/ = or 1 we 

define 

n-l 

fa(k)^l[^aX^'k). (4) 
/=0 

It is clear that fa{k) — fa{—k) and we can check that 

J2fa{k)f;{k)^Sa,b. (5) 

k 

A more detailed analysis shows that fa{k) is non-zero only when k — ±k[a] for a 
function k[a\ which will be specified. More precisely, it is possible to write 

fa{k) = UjSk^k[a] + '^*Sk-k[a] ■ (6) 

This relation defines a one-to-one map k[a]. In the following, a will also be treated 
as an integer modulo 2". We can now expand 

2" 

= J]c«/«(fc) . (7) 
a=l 

With this and Eqs. (5) and (6) we can rewrite 

a=l 

where g{a) = g{k[a]). By construction, the Ca are real field variables. For conve- 
nience we shall also use their complex form 

1 i 

<ya = ^{Ca + C-a) + -{Ca - C-a) ■ (9) 

We can now explain the correspondence between this relabehng and Wilson's 
cell decomposition. Clearly, if = 1) fa{k) is supported in the high momen- 
tum region. More precisely, the 0-th shell, i.c, the one integrated first in the 
RG procedure, consists in configurations which can be expanded in terms of the 
/i+ai2+... (^)- Similarly, the modes corresponding tho the l-th shell are made out 
of the the /2i+a,+i2'+i+....(^)- 
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3. The Hierarchical Approximation 
and its Systematic Improvement 



In the previous section, we have introduced new field variables Ca corresponding 
to the l-ih shell when a can be divided by 2^ but not by 2'+^. This property is 
not affected if a is multiplied (modulo 2") by any odd number. Note also that 
the odd numbers form an abelian group with respect to the multiplication modulo 
2". The orbit of this group within the integers modulo 2"^ are precisely the sets of 
numbers that we have put in correspondence with the shells. The representations 
of this group have been studied and classified. The results can be found in a book 
by Taibleson!^' 

In order to use Taibleson results, we have to embed the labels introduced 
above and denoted a, in the 2-adic integers. Such a technique has already been 
used, for instance, to discuss the random walk representation of the hierarchical 
model!*' When a can be divided by 2^ but not by 2''^\ we say that the 2-adic 
norm, noted |a|2, is In the infinite volume limit, or in other words when n 
tends to infinity, the multiplicative group of the odd numbers is called the 2-adic 
units. The representations of this group will be denoted Ilg. This means that if ui 
and U2 are 2-adic units, then Yls{uiU2) — Tls{ui)^s{u2) ■ The label s specifies the 
representation in a way which will be explained below. 

It is easy to construct explicitly the representations II^. A 2-adic unit can 
be written '^'as u = ±Exp{4:z) where z is a (2-adic) integer and Exp the 2-adic 
exponential. Ils{u) is even or odd under multiplication by —1. On the other hand, 
z is an additive parametrization and the z dependence of Us will be of the form 
where q is an odd integer and r a positive integer. Taibleson calls r + 2 the 
degree of ramification. In summary, the label s is a short notation for the parity, 
r and q an odd integer modulo 2^. 

We can now use these representations to expand the the kinetic term function 
g{a) in each shell. For a given shell I, the a have the form 2'ii (so \a\2 — 2~^) and 
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we can write, 



1 1 S 



(10) 



s 



At finite volume, i.e., at finite n, the units are understood modulo 2""' and con- 
sequently the sum over the representations s is restricted to r < n — Z — 2. The 
numerical coefficients g^^ are easily calculable using the orthogonality relations 
among the representation. 

The hierarchical approximation is obtained by retaining only the trivial rep- 
resentation in the expansion (10). In this approximation, and using the definition 
introduced in Eq.(9), the action reads 



After a Fourier transform, we obtain a hierarchical model having the general form 
written by Dyson in Ref.[3], as discussed below. 

The classification of the representations of the 2-adic units suggest that we 
improve the hierarchical approximation by taking into account the additional terms 
in the expansion (10) order by order in the degree of ramification. Intuitively, 
this corresponds to the fact that the degree of ramification measures the "power of 
resolution" of the representation. Numerically, this works quite well. For instance, 
for g{k) — 1 — cos{^^), we obtain with n — 4 and I — 0, g^^ — 1.63, g_Q — 
0.31, g_^i — 0.25 and g_i — —0.05. The subscripts denote the parity and r 
respectively. Obvious indices have been skipped. 

The actions (8) and (11) take a more familiar form after Fourier transformation 
in a. In general, the l-th momentum shell is responsible for interactions among the 
averages of the Fourier transformed fields inside boxes of size 2^ In the hierarchical 
approximation, the interactions depend only on the position of these boxes inside 
boxes of size 2^+^. When the corrections are introduced up to r = rmax, the inter- 
actions depend on the position inside boxes of size 2'+^+^'"'"^ . If local interactions 




a 
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in these new fields are introduced, then it is easy to generahze the numerical or 
diagrammatic '^"'treatment for small rmax- 

4. The Improvement of the Hierarchical 
Approximation As a Symmetry Breaking Problem 

It is important to realize that the hierarchical approximation of S given in 
Eq.(ll) is invariant under the transformation 

C^a CTua (12) 

for any odd number u. When the other terms of the expansion are incorporated, 
this symmetry is broken by each term in a definite way. This allows us to use Ward 
identities techniques. We discuss the simplest case below. 

Suppose we want to calculate the 2-point function, using the perturbative 
expansion described in the previous section. First we use the new variables and 
the inverse of the map k[a\ defined in section 2 to write 

i^k^-k) = {cl^k]) ■ (13) 

In the hierarchical approximation, the value of this expression depends only on the 
momentum shell specified by \a\2- In other words, 

{cla)o - (Ca>o (14) 

where ( )q means that the quantity is evaluated, at order zero, or in other words, 
in the hierarchical approximation. 

Suppose that we now include a correction Sg{a) to the approximation of g{a). 
Then in first order in this perturbation, we recover the a momentum dependence 
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within the shells given by 

(cla), = (cl), - ^ J2(^~9{ub) - Sm) (clcl), (15) 

b 

In the model considered here, the (cjC^)q contribution can be evaluated straight- 
forwardly and wc can check that wc recover the first term in the expansion of 
{l/g{ua)) — {l/g{a))). The important point is that the corrections are evaluated 
using the unperturbed action. It is clear that similar methods can be used for 
higher point functions and in the interacting case. 

5. Conclusions and Perspectives 

We have proposed a systematic improvement of the hierarchical approximation 
for the gaussian model in one dimension. This procedure exploits in an optimal way 
the symmetries of the approximation and make sense as a perturbative expansion. 

An extension to D-dimension can be constructed easily by noticing the ap- 
proximate correspondence between the integration over successive shells and the 
block-spin method. On a D-dimensional cubic lattice, we can decompose the block- 
spin procedure into D steps (one in each directions). This suggests the introduction 
of functions 

■^i{ki,...kD) =-^o{kj) (16) 
and the replacement of the definition (4) by 

faih kD) = K,{ki) <_,(M^L(2^i) <.-.(2M (17) 

This procedure keeps the computational advantages of the D — 1 case. 

Our general objective is the construction of a perturbative expansion for which 
the RG method and the e-expansion are very well understood in the 0-th order 
approximation. A first test of the reliability of the method will be the calculation 
of the first corrections to the critical exponents. In the longer term, we also expect 
the method to shed some light on Polyakov's conjecture for the 3D Ising model!"' 
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